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I Abstract: In supersymmetric models with gravity- mediated SUSY breaking, universality 

of soft SUSY breaking sfermion masses mo is motivated by the need to suppress unwanted 
flavor changing processes. The same motivation, however, does not apply to soft break- 
^ ■ ing Higgs masses, which may in general have independent masses from matter scalars at 

the GUT scale. We explore phenomenological implications of both the one-parameter and 
two-parameter non-universal Higgs mass models (NUHM1 and NUHM2), and examine the 
parameter ranges compatible with QcDAth 2 , BF(b — > 57) and (g — 2) M constraints. In 
contrast to the mSUGRA model, in both NUHM1 and NUHM2 models, the dark matter 
A-annihilation funnel can be reached at low values of tan /?, while the higgsino dark matter 
annihilation regions can be reached for low values of mo. We show that there may be 
observable rates for indirect and direct detection of neutralino cold dark matter in phe- 
nomenologically aceptable ranges of parameter space. We also examine implications of the 
NUHM models for the Fermilab Tevatron, the CERN LHC and a = 0.5 - 1 TeV e+e" 
linear collider. Novel possibilities include: very light ur, cr squark and S-l slepton masses 
as well as light charginos and neutralinos and H, A and Higgs bosons. 
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Supersymmetric Standard Model. 
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1. Introduction 

The minimal supergravity (mSUGRA) model |l[ provides a convenient and popular tem- 
plate for exploration of many of the phenomenological consequences of weak scale super- 
symmetry [01. In mSUGRA, it is assumed that supersymmetry is broken in a hidden 
sector of the model, with SUSY breaking communicated to the visible sector via gravi- 
tational interactions. The qualifier "minimal" in mSUGRA refers to the assumption of a 
flat Kahler metric, which leads to universal tree level scalar masses at some high energy 
scale, usually taken to be Q = Mqvt- The universality assumption ensures the super-GIM 
mechanism [||, which suppresses unwanted flavor-changing neutral current effects. An at- 
tractive feature of this framework is that electroweak symmetry can be radiatively broken 
(REWSB). This allows one to eliminate the superpotential parameter in favor of M%, 
and the low energy phenomenology is then determined by the well-known parameter space 

mSUGRA : mo, m^, Aq, tan/3, and sign(fi). (1-1) 
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Here mo is the common GUT scale scalar mass, is the common GUT scale gaugino 

mass, Aq is the common GUT scale trilinear term, tan (5 is the weak scale ratio of Higgs 
field vacuum expectation values, and \i is the superpotential Higgs mass term. We take 
m t = 178 GeV throughout this paper. 

The mSUGRA model has been criticized because the assumption of universal scalar 
masses is ad hoc and does not follow from any known symmetry principle Q. While it is 
possible to invoke an additional global U (N) symmetry for the (gravitational) interactions 
of the N chiral supermultiplets, this symmetry is clearly not respected by superpotential 
Yukawa couplings, and radiative corrections involving these Yukawa interactions can lead 
to large deviations from the universality hypothesis ||. 

The assumption of equality of scalar masses receives partial support in Grand Unified 
Theories (GUTs). For instance, in 50(10) SUSY GUT models, all matter superfields of 
a single generation belong to a 16 dimensional spinor representation tp(16) of 50(10), 
and their mass degeneracy is guaranteed if SUSY breaking masses are acquired above the 
5*0(10) breaking scale. If the mechanism by which matter scalars acquire SUSY breaking 
masses is generation blind, universality of matter scalar masses would then obtain. In the 
case of minimal 50(10) SUSY GUTs, the two MSSM Higgs doublet superfields H u and H d 
belong to the same 10 dimensional fundamental representation ^(10), so the corresponding 
SUSY breaking scalar mass terms would not be expected to be the same as those of the 
matter scalars. The phenomenologically desirable super-GIM mechanism would be ensured 
by requiring a U{3) symmetry amongst the different generations. In practice, the amount 
of degeneracy needed is greatest for the first two generations where FCNC constraints 
are the strongest, while the corresponding constraints for the third generation are rather 
mild ||. The need for generational degeneracy can be further reduced if one invokes as well 
a degree of alignment between squark and quark mass matrices, or a (partial) decoupling 
solution to the SUSY flavor problem. 

In this paper, we will maintain degeneracy amongst matter scalars at scales Q ~ 
Mgut, but will allow non-universality to enter the model via soft SUSY breaking masses 
for the Higgs scalars. In our analysis, we will differentiate between two cases for the 
non-universal Higgs mass (NUHM) models. Inspired by GUT models where both MSSM 
Higgs doublets are contained in a single superfield, we will first examine the NUHM model 
where rn 2 H ^ = rrijj ^ m,Q at Q = Mqut- In this case, we define the new parameter 

77i0 = sign(m 2 Hu d ) ■ ^J\m 2 H J at the GUT scale. Thus, the parameter space of this one 
parameter extension of the mSUGRA model is given by, 

NUHM1 : mo, 7770, mi/2-, Aq, tan/3 and sign(fj,). (1-2) 

The second case is inspired by GUT models where H u and H d belong to different 
multiplets. The parameter space for this second case is then given by 

NUHM2 : mo, m^, m H d i m i/2i A)> tan/3 and sign(fi). (1-3) 

The conditions of electroweak gauge symmetry breaking allows one to trade the GUT scale 
masses m 2 H and m 2 H for the weak scale parameters [i and tua- 
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We remark that regardless of any theoretical motivation, if any small departure from 
a well-motivated framework such as mSUGRA causes significant differences in the phe- 
nomenological outcome, the new framework is worthy of examination. We will see below 
that enlarging the model parameter space to split off the GUT scale Higgs boson mass 
parameters from those of other scalars leads to significant departures from mSUGRA 
expectations upon the incorporation of the WMAP constraint on the relic density of 
cold dark matter. Motivated by this, our goal here is to explore in detail the phe- 
nomenological consequences of the NUHM1 and NUHM2 models. Before doing so, we 
note that the mSUGRA model has recently been tightly constrained by several measure- 
ments [0| . These include 1.) the combined measurement || of the branching fraction 
BF(b — > S7) = (3.25 ± 0.37) x 10 -4 , 2.) the measurement || of the deviation of the muon 
anomalous mag netic moment Aa M = A(g - 2)J2 = (27 ± 10) x 10~ 10 from the SM predic- 



tion [10], and 3.) the WMAP determination of the relic density of cold dark matter (CDM) 
in the universe [O: McDMh 2 = 0.113 ±0.009. In addition, we invoke the usual constraint 



from LEP2 that charginos should have mass > 103.5 GeV. We remark that there 

could be significant theoretical uncertainties in the evaluation of both Aa^ and, especially 
for large values of tan (3, also BF(b — > 57) so that any inferences from them should be 
interpreted with care. The first two constraints favor mSUGRA models with [i > 0. The 
WMAP constraint restricts the mSUGRA parameter space to lie in one of the following 
regions 0, H, MM- 

• the bulk region at low mo and low m 1( / 2 , where neutralino annihilation in the early 
universe occurs predominantly via ^-channel slepton exchange (this region is now 
essentially excluded by the combination of WMAP QcDAfh 2 bound and the LEP2 
bounds on and m/j, save where it overlaps with the stau co-annihilation region), 

• the stau co-annihilation at low mo but almost any mi/2 value, where — m^ ]i~6| ], 
or the stop co-annihilation region for special values of Aq where ~ m^ |17]] , 

• the hyperbolic branch/focus point region (HB/FP) at large mo, where becomes 



small, and the neutralino develops a significant higgsino component [18, IS, 20 1, and 



the ^4-annihilation funnel at large tan/3, where 2m^ ~ m^, and neutralino annihila- 
tion in the early universe occurs via the broad A and H Higgs boson resonances plfl . 
A light Higgs resonance annihilation region may also be possible at low m 1 / 2 values 
where 2m^ ~ [22|. 



Throughout this work, we use Isajet 7.72 to generate sparticle mass spectra [23|, IsaReD [13 
for the relic density calculation, and the DarkSUSY package [24| for the computation of dark 
matter detection rates. 

At this point, we would like to take note of a variety of earlier studies of models with 
non- universal Higgs masses. SUGRA models with non- universal Higgs masses were first 
studied by Berezinski et al, who focussed on direct detection of neutralino dark mat- 



ter [27] and indirect detection via neutrino telescopes! 28]. Around the same time, direct 



detection of neutralino dark matter in NUHM models was also investigated by Arnowitt 
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and Nath [29|, and subsequently by Bottino et a^.|pO|- These latter papers explored only 
the case of positive squared Higgs masses. Bottino et al. explored direct DM detection 
for cases including negative squared Higgs masses in Ref. Ellis et al. [32] made a 

thorough exploration of parameter space of the NUHM2 model using the fi and rriA vari- 
ables, and investigated direct detection rates in Ref. |33|. Indirect detection via neutrinos 
was investigated by Barger et al. [34] for models with universality and non-universality. 
Recently, both direct |35|] and indirect 36] detection of neutralino dark matter has been 
investigated by Munoz et al. in the context of models with both scalar and gaugino mass 
non- universality. 

Our present study goes beyond these previous works in several respects: 

1) we investigate the more constrained NUHM1 model in Sec. E|, and show for the first 
time that in this minimal (one parametric extension) of mSUGRA model there are always 
two solutions for low relic density: one is neutralino annihilation via heavy Higgs resonance 
even at low values of tan while the other is neutralino annihilation via higgsino compo- 
nents at low values of mo ; 

2) we find new allowed regions of the NUHM2 model - the light squark/slepton co- 
annihilation regions as discussed in Sec. ^ 

3) we investigate direct and indirect detection of neutralino dark matter, including anti- 
matter, neutrino and gamma ray indirect searches for the new parameter regions mentioned 
above; 

4) for the first time, we consider the implications of the NUHM1 and NUHM2 models for 
collider searches at the Fermilab Tevatron, CERN LHC and ILC linear e + e~ colliders, and 
show how these correlate with the direct and indirect dark matter searches. In this con- 
nection, we also emphasize the sensitivity of the implications of the WMAP measurement 
of QcDMh 2 for collider expectations to the underlying framework. In particular, we show 
that inferences valid in the mSUGRA model may simply be invalid in the extended NUHM 
framework. 

The remainder of this paper is organized as follows. In Sec. |2|, we explore the allowed 
regions of the NUHM1 model which was first studied in Ref. | |25| ]. We will find that even 
for low values of mo, raising the ratio m^/mo brings us into the low region where the 
relic density is in accord with the WMAP allowed range; this is quite unlike the situation 
in mSUGRA where the higgsino annihilation region occurs only at multi-TeV values of 
mo- I n addition, lowering the ratio m^/mo into the range of negative values decreases the 
value of rriA until the ^4-annihilation funnel is reached. In this case, the ^4-funnel region 
can occur at any tan (3 value where an acceptable spectrum can be generated. 

We introduce and outline in Sec 2.3 the computation of direct and indirect dark matter 
rates. We make use of consistent halo models in the attempt to systematically compare 
the reach in all various detection channels. We find enhanced signal rates for direct and 
indirect detection of neutralino cold dark matter in these WMAP-allowed regions [26]. 

In Sec. 2J5, we explore some unique consequences of the NUHM1 model for collider 
searches. In the higgsino region of the NUHM1 model, charginos and neutralinos all become 
rather light, and more easily accessible to collider searches. In addition, lengthy gluino 
and squark cascade decays to the various charginos and neutralinos occur, leading to the 
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possibility of spectacular events at the CERN LHC. In the A- funnel region, the A, H and 
Higgs bosons may be kinematically accessible to searches at the International Linear 
Collider (ILC) or at the CERN LHC, and may also be present in gluino and squark cascade 
decays. 

In Sec. |3|, we explore the NUHM2 model. In this case, since [i and tua can now be 
used as input parameters, it is always possible to choose values such that one lies either 
in the higgsino annihilation region or in the ^-funnel region, for any value of tan (3, mo or 
mi/2 that gives rise to a calculable SUSY mass spectrum. In the low \i region, charginos 
and neutralinos are again likely to be light, and possibly accessible to Fermilab Tevatron, 
CERN LHC and ILC searches. If instead one is in the A-annihilation funnel, then the 
heavier Higgs scalars may be light enough to be produced at observable rates at hadron or 
lepton colliders. In addition, new regions are found where consistency with WMAP data 
is obtained because either ur, cr squarks or left- sleptons become very light. The ur and 
cr co-annihilation region leads to large rates for direct and indirect detection of neutralino 
dark matter, and is in fact already constrained by searches from CDMS2. We present a 
summary and our conclusions in Sec. ||. 



2. NUHM1 model 
2.1 Overview 

In this section, we investigate the phenomenology of the NUHM1 model, wherein the Higgs 
masses m 2 H ^ = fn 2 H = sign(m^) • |m^| 2 at Q = Mqut, with m 2 ^ tuq. We first note that 
the parameter range for m^ need not be limited to positive values at the GUT scale, and 
that, indeed, to achieve radiative EWSB, m 2 ^ must evolve to negative values. Indeed, 
negative squared Higgs mass parameters (at Q = Mqut) are predicted in the SU(5) fixed 



point scenario of Ref. [37 |. 



We also note that some authors impose so-called GUT stability (GS) bounds 

m 



■■? Hu {M GUT ) + ^{M GUT ) > and (2.1) 
m 2 Hd {M GUT ) + fi 2 {M GUT ) > , (2.2) 



to avoid EWSB at too high a scale. The reliability of these bounds is debated in Refs. 
here we will merely remark on regions of parameter space where they occur, and leave it 
to the reader to decide whether or not to impose them. 

We show one of the critical aspects of the NUHM1 model in Fig. [l], where we plot the 
values of //, tua, tu^ and 2m^ versus m^ while fixing mo = mi/2 = 300 GeV, with Aq = 0, 
tan/3 = 10 and \i > 0. The region to the left of the dot-dashed line indicates where the GS 
bound fails. The curves terminate because electroweak symmetry breaking is not obtained 
as marked on the figure: in fact, on the right, where becomes small, the chargino mass 
falls below the LEP bound just before the EWSB constraint kicks in. The black curve 
denotes the value of the fi parameter, which takes a value of [i = 409 GeV for ms = 300 
GeV (the mSUGRA case). The parameter [i becomes much larger for m^ < —mo, and 
much smaller for m^ > tuq [5^]. The region of small \x is of particular interest since in 
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that case the lightest neutralino develops substantial higgsino components, and leads to a 
relic density which can be in accord with the WMAP determination. In contrast, in the 
mSUGRA model the higgsino-LSP region occurs in the HB/FP region, which occurs at 
very large m,Q values of order several TeV (depending somewhat on the assumed value of 
mt). The HB/FP region has been criticized in the literature in that the large mo values 
may lead to large fine-tunings |^] (for an alternative point of view, see Refs. [18] and [|lS|]), 
This illustrates an important virtue of the NUHM1 model: the higgsino annihilation region 
may be reached even with arbitrarily low values of mo and mj^, provided of course that 
sparticle search bounds are respected. 

We also see from Fig. [l] that the value of rriA can range beyond its mSUGRA value for 
large values of m^, to quite small values when becomes less than zero. In particular, 
when rriA ~ 2m^ , neutralinos in the early universe may annihilate efficiently through the 
A and H Higgs resonances, so that again QcDMh 2 may be brought into accord with the 
WMAP result. In the mSUGRA model, the ^-annihilation funnel occurs only at large 
tan (5 ~ 45 — 55. However, in the NUHM1 model, the A-funnel region may be reached even 
for low tan/? values if is taken to be sufficiently negative. 
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Figure 1: Plot of fi, ttia and vs. for m — 300 GeV, m 1 / 2 = 300 GeV, A = 0, tan/3 = 10 
and m t = 178 GeV. We take \i > 0. 



To understand the behavior of the /i parameter and ttia in the NUHM1 model, we first 
examine the renormalization group equations (RGEs) for the soft SUSY breaking Higgs 
squared mass parameters. Neglecting the Yukawa couplings for the first two generations, 
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these read: 

^ = is? H 9 ? M > - 3 ^ + + 3/ « x *) • < 2 - 3) 
= il? H 9?M? - % ' M| - + 3/ » Xb + • (24) 

where t = log(Q), /t,6,r are the t, b and r Yukawa couplings, and 

X t = m 2 Qa + m? fl + m 2 Hu + A?, (2.5) 

X 6 = m 2 Q:j + + m 2 Hd + ,4 fe 2 , (2.6) 

V T = m| 3 + m~ fl + m 2 Hd + and (2.7) 

S = m 2 Hu - m 2 Hd + Tr [m 2 Q - m| - 2m 2 , + m 2 D + m|] . (2.8) 

The term S 1 is identically zero in the NUHM1 model, but can be non-zero in the NUHM2 
model. 

For small-to- moderate values of tan/3, / ( > /j, f T , and so the RGE terms including 
Xt usually dominate the Xf, and X T terms. The RGE terms including Yukawa couplings 
occur with overall positive signs, which results in driving the corresponding soft Higgs 
boson mass squared parameters to smaller (and ultimately negative) values at the low 
scale. Indeed, this is the familiar REWSB mechanism. Since Xt 3 m2 H u -> a large, positive 
value of rritf) > tuq results in a stronger push of to negative values (relative to that in 
mSUGRA) during the running from Mqut to M wea k, while the evolution of rrijj is rather 
mild. Alternatively, if <C 0, then there exist cancellations within the Xt term which 



results in a milder running of m H from Mqut to M wea k. Indeed, as shown in Ref. 25] 



Am 2 Hu d = m 2 Hu d (NUHMl) - rn% u d (mSUGRA) 

satisfies 

Am^(weak) ~ Am^ (GUT) x e~ Jt , 

where 



J t = 7T2 I dt ft > °> 



3_ 

in 1 

with ft being the top quark Yukawa coupling. We see that Am^ maintains its sign, 
but reduces in magnitude under evolution from the GUT to the weak scale. The same 
argument applies for Arrijj , except that the effect of evolution is much smaller because 
fb,r ^ ft except when tan (3 is very large. The situation is illustrated in Fig. |2| where we 
plot the running of rrijj and q rn\i d from Mqut to M wea k using the same model parameters 
as in Fig. [j], except for three choices of = 500 GeV, 300 GeV (mSUGRA case) and 
—700 GeV. In these cases, the weak scale values of m 2 H are —(251 GeV) 2 , —(407 GeV) 2 
and —(732 GeV) 2 , respectively, while the corresponding weak scale values of m 2 H are 
(527 GeV) 2 , (348 GeV) 2 and -(672 GeV) 2 . 

The tree level minimization condition for EWSB in the MSSM is 

_ m^-m 2 ^ tan 2 /? _ M 2 

~ (tan 2 /?-l) 2 ' ^ iUj 
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m =300GeV, m 1/2 =300GeV, tan$=10, A =0, \i>0, m t =178GeV 
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Figure 2: The evolution of m 2 Hu and m 2 Hd from Q = M GUT to M weak for M$ = 300 GcV (mSUGRA 
case), 500 GeV and -700 GeV. We also take m = 300 GeV, m 1/2 = 300 GeV, A = 0, tan/3 = 10 
and mt = 178 GeV, with /i > 0. 



For moderate to large values of tan/3 (as favored by LEP2 Higgs boson mass constraints), 
and |m#J S> Mz, fJ- 2 ~ —m 2 Hu . Thus, we see that in the case of large negative 
values, we would expect a large parameter, whereas for large positive values, m# 
is barely driven to negative values, and we expect a small parameter. Within the same 
approximation, the tree level pseudoscalar Higgs mass ttia is given by 



m H u + ™H d + - m H d 



(2.11) 



For large negative values of m^, the weak scale values of and vnr H are both negative, 
and can cancel against the 2fi 2 term, yielding small pseudoscalar masses. Meanwhile, for 
large positive values of m^, m 2 Hd ~ signim^rn 1 ^ while m? Hu is small, but negative. In this 
case there is no cancellation in the computation of m\, and thus we expect tua to be large, 
as shown in Fig. [|. 

In Fig. ||, we show how different sparticle masses vary with for the same parameter 
choices as in Fig. |]. Most sparticle masses are relatively invariant to changes in ms. One 
exception occurs for the W\ and Z\% masses, which become small when /j, ~ M2, and the Z\ 
becomes increasingly higgsino-like. The other exception occurs for the i\ <i and b\^. masses. 
In this case, the Q3 = &l) and tR running masses also depend on terms including f 2 X t . 
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m =300GeV, m 1/2 =300GeV, tan\i=10, A =0, (j. >0, m t =178GeV 
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Figure 3: Various sparticle masses versus in the NUHM1 model for mo = 300 GeV, m 1 / 2 — 300 
GeV, A = 0, tan (3 = 10 and m t = 178 GeV. We take n > 0. 



Thus, when Xt is small (for < mo), the diagonal entries in the top and bottom squark 
mass squared matrices are not as suppressed due to top Yukawa coupling effects as in the 
mSUGRA case. In contrast, for large positive values of m^, Xt is large and these soft 
masses are more suppressed resulting lighter third generation squarks. These expectations 
are indeed born out in Fig. ||[ 

2.2 NUHM1 model: parameter space 

The mSUGRA parameter space point we have used for illustration so far, 

(m , mi/2, Aq, tan (3 = 300 GeV, 300 GeV, 0,10) 

with [x > and m t = 178 GeV, is excluded since ^^^i 2 = 1.2. However, by extending the 
parameter space to include as in the NUHM1 model, these parameter values are allowed 
for an appropriate choice of m^. As an example, in Fig. || we use mo,mu2 and Aq as in the 
mSUGRA parameter set above, and scan over tan/3 and m^ values. In plotting points, 
we construct a x 2 value out of the three quantities Q^/i 2 , BF(b — > S7) and Aa^, where 



X 



El 



1 2 — , where Xi is the predicted value, \ii is the measured value, and cij is the 

error on the ith measured quantity. In constructing the x 2 value, we only use the WMAP 
upper bound (thus, points with f^/t 2 < 0.113 do not contribute to the \ 2 ) to allow for the 
possibility of mixed cold dark matter, where for instance a portion of dark matter might 
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m D =300GeV, m 1/2 =300GeV, A =O f \i>Q, m t =178GeV 
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Figure 4: Ranges of \/x^ in a scan over tan f3 and values in the NUHM1 model for fixed 
toq = mi/2 — 300 GeV, A n = 0, > and m t = 178 GeV. The green region corresponds to low 



values of yx ~ V3, while the red region has VX^ 
intermediate values. 



5, with the yellow region corresponding to 



consist of, say, axions 1 . Green points have low yx ~ "V^ values, and agree well with the 
central values of each of these measurements. Red points have large -\/x 2 ~ ~ 5, and are 
excluded. Yellow points have intermediate values of \fx*. We show six frames illustrating 
various correlations amongst parameters. In frame a) showing vs. tan /?, we see the 
green/yellow A-annihilation funnel for ~ —0.8 TeV, which occurs for every tan (3 value. 
We also see at ~ 0.6 TeV the appearance of the higgsino region, corresponding to the 
HB/FP region of the mSUGRA model. While the relic density is in accord with WMAP 
in this region, as tan/3 increases, BF(b — > s'y) and Aa^ also increase, so that the higgsino 
region becomes increasingly disfavored for large tan/3. Frame b) shows the tan/3 vs. \x 
correlation, where we see that indeed /x is small in the higgsino region, and large in the 
A-funnel. Frame c) shows the vs. tha correlation, and indeed we see large values of 
rriA in the higgsino region, while itia ~ 250 GeV in the A-funnel. The remaining three 
frames show (2m^ — m a) I'm A vs. tern (3, and \i respectively, which explicitly displays 
the ^4-annihilation funnel against the input parameters and \x. 

As we have already mentioned, there is still some debate on the range of the SM 
prediction for gu, as well as the MSSM prediction for BF(b — > 57). As a result, the x 2 

1 A low thermal relic abundance may also be compatible with a fully supersymmetric dark matter scenario 
provided non-thermal production of neutralinos or cosmological enhancements of the thermal relic density 



occur: see Sec 2.3 
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Figure 5: Ranges of fi/i 2 , 5<x M and BF{b — > S7) for the NUHM1 models scanned in Fig. 
parameter point falls outside the ranges shown, it is not plotted in this figure. 



|. If a 



values in Fig. |] should be interpreted with some care. To facilitate this, we show the ranges 
of 17/i 2 , 5a^ and BF(b — > S7) for the same set of NUHM1 models as in Fig. || If a parameter 
set in Fig. |I| yields a value of these quantities that is outside the range shown, then the 
point is not plotted. From this plot we see why only the low tan (3 portion of the higgsino 
region at low /i gives low x 2 m Fig. |j: at high tan /3, the BF(b — > 57) is quite high, while the 
value of Aa^ is quite low. We see, however, a large swath of yellow with tlcDMh 2 < 0.094 
at large values of tan f3 in the upper frames of Fig. ||. This occurs because the A and 
H bosons become light and relatively wide, leading to a resonant enhancement (even off 
resonance) in the neutralino annihilation cross section and a corresponding reduction in 
the relic density. In this case, there must then be some other new physics that brings the 
CDM density up to the WMAP value 2 . 



2 We note that it is possible that this new physics associated with the non-LSP components of dark 
matter may also yield (possibly non-calculable) new contributions to both 8a M and BF(b — » S7), so that 
it may be premature to unequivocally exclude this parameter space region at large values of tan /3 because 
these quantities are not in agreement with their measured values. 
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Figure 6: The a/x^ value in the m vs. m\n plane for A — 0, tan/5 = 10, fi > and m t = 
178 GeV. In frames a), 6), c), and <f) we take m^/mo = 1, —2.5, 1.1 and 1.5, respectively. The 
blue region is excluded by LEP2. 

In Fig. [| we show the mo vs. m^ 2 parameter space plane for the NUHM1 model, 
with Aq = 0, tan/? = 10, \i > 0, m t = 178 GeV and a) = (mSUGRA case), 6), 

= — 2.5mo, = l.lmo and = 1.5mo- The black regions are excluded by lack of 
REWSB (right hand side) and because the stau is the LSP on the left hand side. The blue 
shaded region is excluded by the LEP2 constraint that < 103.5 GeV. The remaining 
parameter space is color coded according to the s/x* value, and indeed we see that most 
of parameter space is excluded. The mSUGRA case of frame a) shows the HB/FP region 
at mo ~ 8 — 10 TeV, while the stau co-annihilation is squeezed against the left edge of 
the allowed parameter space. In frame b) for a large negative value of ms, we see that a 
narrow allowed region now cuts through the middle of the parameter plane. This is the 
yl-annihilation funnel, which is much narrower than in the mSUGRA case at large tan (3, 
since now the A-width is relatively small: typically ~ 1 GeV. Note that the range of mo 
extends only to 2 TeV, since for larger values of mo, rn 2 A < and REWSB is violated. Since 
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\x remains large, there is no higgsino LSP region along the right-hand edge of parameter 
space. In frame c) for = l.lmo, we see that the parameter ranges only to 3 TeV, 
since now the right-hand side is excluded by u? < 0. This leads to a higgsino LSP region 
which is shaded yellow, which begins at ttiq ~ 1 TeV for low mi/2, an d explicitly shows 
that the higgsino LSP region can occur even for relatively light scalar masses. Frame d) 
for ms = 1.5mo shows that the higgsino region has moved to even lower mo values, which 
are below 3 TeV even for mi/2 as high as 2 TeV. 



NUHM1: tan$=35, m o = -2.5111^ n >0, A Q =0, m t =178 GeV 
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Figure 7: The a/x^ value in the mo vs. m 1 / 2 plane for Aq = 0, tan/3 = 35, /i > and m t — 178 
GeV, with m</, = — 2.5toq- The blue region is excluded by LEP2. 



In Fig. 0, we show the \px 1 ~ value in the tuq vs. m^j plane for Aq = 0, tan/3 = 35, 
fi > 0, = — 2.5mo and mt = 178 GeV. This plane includes the region of lowest x 2 value 
as indicated in Fig. ^j. Here we see a well-defined ^4-annihilation funnel for the case of 
tan /3 = 35, where the lower portion gives excellent agreement to the measured WMAP relic 
density, the branching fraction BF(b — > S7) and the muon anomalous magnetic moment 
Aa^. Excellent agreement is also obtained in this case for the stau co-annihilation region 
for mi/2 ~ 350 - 600 GeV. 
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2.3 Dark matter detection: overview and methodology 

Two major issues enter the evaluation of the prospects for neutralino dark matter detection 
and of assessing the relative effectiveness of the various direct and indirect techniques that 
have been proposed: 

1. Since the detection rates critically depend on the dark matter halo profile of our own 
galaxy, and a wealth of observational data and theoretical constraints are available, 
the halo models one resorts to must be consistent with all information we have about 
the Milky Way. Moreover, since different detection techniques rely on the local dark 
matter density distribution and on the velocity distribution, a fair comparison among 
them may be carried out only provided the two distributions are self- consistently 
computed, not only locally, but throughout the whole halo 3 . 

2. A comparison of the various techniques must rely on quantities which provide infor- 
mation on the relative strength of the expected signal with respect to the projected 
experimental sensitivity, taking into account the background. Toward this end, we 



use what have been dubbed Visibility Ratios (VR) |44| , Eq , 46, 25], which will be 
defined below for each experiment. A VR is simply a signal-to-sensitivity ratio: when 
VR>1, the signal calculated using a specified model is expected to be detectable over 
backgrounds with the particular experimental setup; in case VR<1, the signal will 
lie below the sensitivity of the considered detection facility. The locus of points at 
VR=1 outlines the actual reach contour of any particular detection technique. The 
relative magnitudes of the VR's for different dark matter search experiments provides 
a direct comparison between themp^j. 



As far as item 1. is concerned, we follow here the strategy outlined in Ref. |44[] (the 
reader is referred to Ref. [It], |45|] for more details). Motivated by N-body simulations 
of hierarchical clustering in dark matter cosmologies, the starting point is to take into 
account the correlation between the virial mass of a given galaxy M v i r and its concentration 
parameter c v i r J5C| ], and then to provide a description of the dynamics of the baryon infall 
and of its back-reaction onto the dark halo. Two extreme regimes have been outlined, the 
first one being motivated by a large angular momentum transfer between dark matter and 
baryons, while the second one by supposing that baryons settled in with no net transfer of 
angular momentum to the dark component. 

In the first scenario, the central cusp in the dark matter halo, as seen in numerical 



simulations, is smoothed out by a significant heating of the cold particles [51], leading to 



a cored density distribution, which has been modeled by the so called Burkert profile [|52j , 

PbW = n I / Vihrt i w - (2 - 12) 

(1 + r/a) (1 + [r/aY) 



3 We recall that, while the gamma ray flux from the galactic center only depends on the dark matter 
density in a very small, and poorly known region, the antimatter rates vary with the dark matter density 
in a more extended portion of the halo; finally the flux of neutrinos from neutralino annihilations in the 
sun and the direct detection rates essentially depend only on the density and velocity distributions of the 
dark halo within the solar system. 
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Here, the length scale parameter has been set to a = 11.7 kpc, while the normalization 
p° B is adjusted to reproduce the local halo density at the Earth position to PB(ro) = 
0.34 GeV cm -3 f48|| . We refer to this model as to the Burkert Halo Model. It has been 
successfully tested against a large sample of rotation curves of spiral galaxies [|53|]. 

In the second scenario that we consider, baryon infall causes a progressive deepening 
of the gravitational potential well at the center of the galaxy, resulting in increasingly 
higher concentration of dark matter particles. In the circular orbit approximation [54, 



55], this adiabatic contraction limit has been worked out starting from the N03 profile 



proposed in Ref. [56]; the resulting spherical profile, which has no closed analytical form, 
roughly follows, in the inner galactic regions, the behavior of the profile proposed by Moore 
et oZ.,|57], approximately scaling as r -L5 in the innermost regions, and features a local 



dark matter density pno^{tq) = 0.38 GeV cm 3 . We dub this setup as the Adiabatically 
Contracted N03 Halo Model. 

The parameters for both models have been chosen to reproduce a variety of dynamical 
information, ranging from the constraints stemming from the motion of stars in the sun's 
neighborhood, total mass estimates from the motion of the outer satellites, and consistency 
with the Milky Way rotation curve and measures of the optical depth toward the galactic 



bulge [48, Finally, the local velocity distributions of the two spherical profiles have 
been self-consistently computed using the formalism outlined in Ref. ]^8| , thus allowing for 
a reliable comparison between direct and indirect detection techniques. Both models have 
been included in the latest public release of the DarkSUSY package |24]] . 

Turning to item 2., we begin by defining the Visibility Ratios (VR) mentioned above 
for each dark matter detection technique. For direct dark matter detection, in view of 
the fact that spin-dependent searches have been shown to be sensitive to scattering cross 



sections typically much larger than those predicted within the MSSM [45], we focus here 



only on spin-independent (SI) searches. The relevant quantity is the neutralino-proton 
SI scattering cross section, er~ } . The VR is simply defined as the ratio between the 

expected cr- 1 from a given supersymmetric model and the corresponding experimental 

Zip 

sensitivity cTg* p at the given neutralino mass, for the particular dark matter halo profile 
under consideration, taking into account the corresponding self-consistently computed local 
velocity distribution. 4 We consider here two benchmark direct detection experiments, 
namely CDMS-II |5^] (stage-2 detectors) and the proposed XENON 1-ton facility (stage-3 
detectors) |6^, 31 1. 



To assess the sensitivity of a km 2 -size detector such as IceCube designed to detect 
high energy muons, we compute the expected muon flux from neutralino annihilations in 
the sun 5 above a threshold of 1 GeV. The low threshold was chosen to allow a comparison 
of the expected muon flux with the existing upper limits from the SuperKamiokande ex- 



4 The sensitivity of direct detection experiments may vary over more than two orders of magnitude for 
the typical range of neutralino masses in the MSSM. It also depends (less strongly) on the halo profile. In 
view of this sensitivity, we must be careful to use the corresponding value of a^ p (m^ ) in the computation 
of the VR. 

5 We also evaluated the muon flux from the Earth, but always found far smaller fluxes than those expected 
from the sun. 
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periment (which indeed has a much lower energy threshold than IceCube) and to ensure 
that the model is not already excluded. The predicted muon flux is then compared with 
the sensitivity of IceCube. Since this detector has a higher energy threshold ~ 50 GeV, 
we corrected for the threshold mismatch using the corresponding projected sensitivities 
worked out in Ref. 63]. We have checked that the projected maximal sensitivity to the 
muon flux obtained using this procedure is in good agreement with that obtained by in- 
tegrating above the actual IceCube threshold. We also incorporated the dependence of 
the sensitivity on the soft (e.g. from neutralino annihilation to bb) and hard (e.g. from 
annihilation to W^ + iy _ , Z°Z° and tt) neutrino spectra — the sensitivity being smaller in 
the case of the softer spectrum. 

Turning to antimatter experiments, we compute the solar-modulated positron, antipro- 
ton and antideuteron fluxes, following the procedure outlined in Ref. [44]. We calculate 
the neutralino annihilation rates to p and n using the Pythia 6.154 Monte Carlo code [64] 
as implemented in DarkSUSY [ 24 1 , and then deduce the D yield using the prescription sug- 
gested in Ref. ||| . The propagation of charged cosmic rays through the galactic magnetic 
fields is worked out through an effective two-dimensional diffusion model in the steady 
state approximation, while solar modulation effects were implemented through the analyt- 
ical force-field approximation of Gleeson and Axford [36|. The solar modulation parameter 
<&F is computed from the proton cosmic-ray fluxes, and assumed to be charge-independent. 
The values of &f we use for antiprotons and positrons refer to a putative average of the 
solar activity over the three years of data taking of the PAMELA experiment 67], which 
will be the first space-based experiment for antimatter searches. 

For antideuterons we consider the reach of the proposed gaseous antiparticle spectrom- 
eter (GAPS) HI proposed to be placed on a satellite orbiting around the Earth, and tuned 
to look for antideuterons in the very low kinetic energy interval from 0.1 to 0.4 GeV |68|, |69|. 
In this energy interval, the estimated background is strongly suppressed ]65|. fr"0| , and un- 
ambiguous evidence of even a single low energy antideuteron could be regarded as a positive 
search result. Since the experiment is expected to be launched close to the maximum of 
the solar cycle, we set the value of the solar modulation parameter &p at the correspond- 
ing value. The resulting sensitivity of GAPS has been determined to be of the level of 
2.6 x 10~~ 9 m~ 2 sr -1 GeV s _1 |3^]. The VR for antideuterons will then be given by the 
integrated D flux over the kinetic energy interval 0.1 < Tjj < 0.4 GeV, divided by the 
GAPS sensitivity. In other words, the VR will give, in this case, the actual number of 
antideuterons expected to be detected by GAPS. 

Regarding antiproton and positron fluxes, it was realized long ago that the low energy 
tails are considerably modified by solar modulation effects and, in the case of antiprotons, 
also by large secondary and tertiary backgrounds [71], so that low energy positrons and 
antiprotons cannot provide a clean test for new physics contributions. However, in view 
of the fact that novel space-based experiments will be able to extend the experimental 
sensitivity to these to the hundreds of GeV range, the best place to look for an antimatter 
signal from neutralino annihilations lies at the high energy end. Moreover, a distinctive 
signature is provided by the clean cutoff of the neutralino-induced antimatter flux corre- 
sponding to energies equal to the neutralino mass, provided this flux is large enough to 
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be detectable. Nevertheless, in general, the energy spectra of antiparticles generally look 
quite featureless: aside from the fact that the antiparticles from neutralino annihilation 
sit on top of a featureless background, solar modulation and propagation effects also tend 
to wash away any features inherent to the spectrum. With some obvious exceptions, the 
antiparticle energy (or rather E/m^ ) where the signal-to-background ratio is largest, is 
sensitive to the composition of the neutralino [72, 47]. Restricting attention to particular 



energy bin(s) may also be misleading because of large uncertainties from secondary and 
tertiary p. A possible exception may be that of neutralinos mainly annihilating into gauge 
bosons, in which case one expects a bump in the positron spectrum at approximately half 



the neutralino mass 73]; this could be especially relevant for neutralinos in the HB/FP 
region. 

To circumvent these problems, we adopt the statistical treatment of the antimatter 



yields introduced in Ref. [44] (an analogous approach has been proposed for cosmic positron 
searches J74j). Motivated by the fact that the signal is much smaller than the background, 
we introduce a quantity which weighs its "statistical significance, summed over the energy 
bins", 

h ~L a ^MEY dE > (2 - 13) 

where (f> s (E) and (fib(E) respectively represent the antimatter differential fluxes from neu- 
tralino annihilations and from the background at antiparticles' kinetic energy E, and 
^min, m ax correspond to the antiparticle's maximal and minimal kinetic energies to which 
a given experiment is sensitive (in the case of the PAMELA experiment @, T m | n = 50 
MeV, T m t x = 270 GeV, T min = 80 MeV and T max = 190 GeV). We compute the primary 
component, (j) s , with the DarkSUSY package, interfaced with a subroutine implementing the 
diffusion and solar modulation models outlined above. The background flux (pi, has been 



calculated with the Galprop package [75], with the same parameter choices employed to 
compute the signal. 

Given an experimental facility with a geometrical factor (acceptance) A and a total 
data-taking time T, it has been shown [[5J] that, in the limit of a large number of energy 
bins and of high precision secondary (i.e. background) flux determination, a SUSY model 
giving a primary antimatter flux (j) s can be discriminated at the 95% C.L. if 

I4<P S )-A-T>( X XT> ( 2 - 14 ) 

where (x 2 )n^° stands for the 95% C.L. x 2 with nb degrees of freedom. For the PAMELA 
experiment, where A = 24.5 cm 2 sr, T=3 years and nb — 60 we get the following discrimi- 
nation condition | |44| ] 

Ws ) > = /Jy. PAMELA, 95% „ 3 2 x 1Q -8 ^-2^.-^-1 (2 15) 

which is approximately valid for both positrons and antiprotons (though in the latter case 
the PAMELA experiment is expected to do slightly better). As a rule of thumb, the 
analogous quantity for AMS-02 should improve at least by one order of magnitude f73fl. 
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We therefore define as VR for antiprotons and for positrons the ratio 



(VR) p ' e = I p f /If' ™ 1V1 ^ A - yo/0 . (2.16) 

Finally, the case for the gamma ray flux from the galactic center is plagued by large 
uncertainties on the very central structure of the Milky Way dark halo. Depending on 
various assumptions on the galactic models and on the physical cut-off in the inner part 
of our galaxy, there might be a spread of various orders of magnitude in the computation 



of the actual gamma-ray nux[|7q, 26 1. This may be written as a product of a purely 
astrophysical quantity describing the propagation of the photons to the detector, and of 
a purely particle-physics quantity (av) 77] describing the source of these photons. The 



former reads, 



where 



<J(0))(Af2) = ^jT J(^)dn, (2.17) 



•'H= tAzZ ■ ( 7^T7^77Z3 ) / P 2 (l>V(r). (2.18) 

line of sight 



.5 kpc V - 3 GeV/cm J 



The attitude we take here is just to extrapolate the halo models we use in the p(r — ► 0) 
limit, and to compute the corresponding (J(0)) for the acceptance Af2 of GLAST. What 
we find is (J(0)> = 7.85 for the Burkert Halo Model, and (J(0)) = 1.55 x 10 5 for the 
Adiabatically Contracted N03 Halo Model. 6 We then compute the integrated gamma 
ray flux above a 1 GeV threshold, (fP , and define the corresponding VR as (ffl/(1.5 x 
1CP 10 cm _2 s~ 1 ), the latter being the corresponding estimated sensitivity of the GLAST 



satellite [|T9 1 . 



2.4 Dark matter detection: the NUHM1 model 



In Fig. H and g we show the VR's for the various experiments that we detailed in Sec. 2.3 
for two representative mSUGRA parameter choices: we fix tan/3 = 10, mo = 300 GeV, 
m 1/2 = 300 GeV in Fig. §, and tan/3 = 20, m = 1000 GeV, m 1/2 = 200 GeV in Fig. g 
(in both cases Aq = 0, sign(/z) > and wi[ = 178 GeV) and show the results as a function 
of m^/mQ. The regions shaded in red are excluded by the LEP2 limits on the mass of 
the pseudoscalar Higgs boson rriA and on the mass of the lightest chargino. The green 
regions indicate, instead, parameter space regions where the neutralino relic abundance 
h? < 0.13, consistently with the WMAP 95% C.L. upper limit on the Cold Dark 
Matter abundance |[(J : agreement with the central value of WMAP is obtained close to 
the boundary of this region. We remind the reader that a VR larger than unity means 
that the signal should be detectable in the particular dark matter detection channel. For 
definiteness, we adopt the conservative Burkert Halo Model. The results we show should be 



6 In the case of the Adiabatically Contracted halo model, this procedure is quite arbitrary, since different 
hypotheses on the dynamics of the central black hole formation might lead to very different predictions for 
the dark matter density in the center of the halo jHf], Q . Because there is an unconstrained extrapolation, 
essentially any flux may be possible. 
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regarded as plausible lower limits, particularly as far as indirect rates are concerned since 
a cuspy inner dark halo would greatly enhance the dark matter detection rates jwl, |7(|. 7 
We see, in Fig. ||, that for all values of vha the signal will be accessible to stage-3 direct 



detection facilities, like XENON 1-ton |60|, while stage-2 detectors (such as CDMS-II^g) 
will be able to probe only the HB/FP region, at large m^. The behavior of the direct 
detection VR's is readily understood. On the one hand, when takes large negative 
values, both tua and the CP-even heavy Higgs boson mass ran decrease; this leads to an 
enhancement of the i-channel H exchange in the neutralino-proton cross section (which 
scales as rnj^, assuming that mA,H 3> On the other hand, when is large and 

positive, the higgsino fraction increases, and so does erf 1 since erf 1 oc (ZhZ Q ) 2 , where 

Zip Zip 

Zh, g respectively denote the higgsino and bino fraction in the lightest neutralino. The 
same behavior for the direct detection VR applies to Fig. |9|; here, once again, stage-3 
detectors will be able to fully explore the parameter space, while only the focus point 
region will be discoverable at stage-2 facilities. 

Turning to the neutralino- annihilations-induced flux of muons from the sun, we see 
that in both figures the resonant annihilation region gives rates which will be various 
orders of magnitude below the projected ultimate sensitivity of IceCube. This is because 
the neutralino capture rate in the sun (which mainly depends on the neutralino-proton spin- 
dependent scattering cross section) is not large enough: in contrast to the spin-independent 
cross section, in fact, the main contribution to the spin-dependent one comes from the Z 
exchange diagram, which is not enhanced by the smaller values of rriH,A- On the other 
hand, a larger higgsino fraction and annihilation cross section yields detectable rates at 
neutrino telescopes for the model considered in Fig. ||, at large m^; for the case addressed 
in Fig. ^ we see that the expected rates only lie less than one order of magnitude below 
the maximal sensitivity of IceCube. This is because for the smaller value of m-i/2) the 
neutralino LSP does not acquire a sufficiently large higgsino component all the way to the 
LEP2 limit. 

As far as other indirect detection techniques are concerned, though we are here con- 
sidering the conservative cored dark matter profile, the large neutralino pair annihilation 
rate {ctv)q in the funnel region yields very large rates in all channels, peaked on the value 
of at which tua — 2m ^ . A considerable enhancement is also seen in the HB /FP re- 
gion. Despite this enhancement, antimatter detection rates might not be large enough to 
be discriminated against the background, especially in the case that we study in Fig. ||. 

We should mention that for ranges of parameters that yield a thermal dark matter 
density smaller than the WMAP central value we do not correspondingly scale down the 
results in Fig. |8| and Fig. |9[ as we would have to for any model where the remainder of the 
dark matter was composed of something other than the SUSY LSP (see e.g. Ref. and 
references therein). In this case, possible additional contributions from these other dark 
matter components would have to be included. We work here, instead, under the hypothesis 
that even within the low thermal relic density parameter ranges, the LSP is all of the dark 
matter, but that there is either additional non-thermal production [6^] , or cosmological relic 



7 This is explicitly illustrated for the NUHM2 Model in the next Section. 
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Figure 8: Dark matter detection Visibility Ratios (i.e. signal-to-projected-sensitivity ratios) for 
various direct and indirect techniques, as a function of the GUT scale non-universal Higgs mass 
parameter to^, for mSUGRA parameters tan [3 = 10, A = 0, sgn// > 0, m = 300 GeV, m 1 / 2 = 300 
GeV, and setting m t = 178 GeV. We take the dark matter distribution to be given by the cored 
Burkert Halo profile described in the text. A Visibility Ratio larger than one means that the signal 
will be detectable over background. 



density enhancement, as envisaged in Ref. [83, S4| for quintessential cosmologies, in Ref. [85] 
for Brans-Dicke-Jordan cosmologies, and in Ref. [ 



87] for anisotropic cosmologies. 
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Figure 9: The same as Fig. || but with mSUGRA parameters tan/? = 20, Aq — 0, u > 0, 
too = 1000 GeV, to 1/2 = 200 GeV. 



2.5 NUHM1 model: Collider searches for SUSY 

While direct and indirect detection techniques that we have just discussed could establish 
the existence of dark matter, collider experiments would be needed to make the link with 



supersymmetry [88]. Collider expectations within the NUHM1 framework can differ from 
corresponding expectations within the well-studied mSUGRA model. In the following dis- 
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cussion, we highlight these differences confining our discussion to NUHM1 model parameter 
sets that satisfy the WMAP bound on fi^/i 2 . 

2.5.1 Fermilab Tevatron 

The most promising avenue for a supersymmetry discovery at the Fermilab Tevatron in the 
case of gravity-mediated SUSY breaking models with gaugino mass unification is by the 
observation of trilepton signals hompp — > W1Z2X followed by W\ — ► lvgZ\ and Z2 — ► 



three body decays, where I = e or |89, 9C|. In the case of the NUHM1 model where m^, is 
taken to have negative values so that neutralinos annihilate via the A and H poles, the only 
effect on the gaugino/higgsino sector is that the magnitude of the \i parameter increases. 
As a result, the HB/FP region of small is absent in frame b) of Fig. so that unlike 
the case of mSUGRA |9(J, probing large values of mo and mj/j via this channel will not be 
possible. Alternatively, if m^ is taken to be large compared to mo so that is small, then 
chargino and neutralino masses can become lighter, which may increase production cross 
sections. Furthermore, the rag — mg mass gap will diminish, which can close "spoiler" 
decay modes such as Z2 — * Z\h, so that the necessary three-body neutralino decays are 
more likely to be in effect, and because tan (5 is not necessarily large, we do not expect 
events with tau leptons to dominate at the expense of e and \i events. Thus, in the large 
m,0 region, we expect improved prospects for clean trilepton signals. Detailed simulation 
would of course be necessary to draw definitive conclusions. 

2.5.2 CERN LHC 

The CERN LHC is expected to begin operation in 2007 with pp collisions at y/s = 14 TeV. 
In most regions of parameter space of gravity-mediated SUSY breaking models, gluino and 
squark pair production is expected to be the dominant source of sparticles at the LHC. 
Since the values of mo and mi/2 determine for the most part the magnitudes of the squark 
and gluino masses, we expect sparticle production rates in the NUHM1 model to be similar 
to those in the mSUGRA model for the same model parameter choices. The reach of the 
CERN LHC in the case of the mSUGRA model has recently been re-evaluated in the 
mo vs. mi/2 plane for various tan j3 values, and assuming 100 fb _1 of integrated luminosity 
in Ref. [MJ. We display this reach contour on frame a) of Fig. 1C where as in Fig. 



a 



we take Aq = 0, tan/3 = 10 and = too (mSUGRA case). We also show the WMAP 
allowed region (% h 2 < 0.13) as the one shaded in green. The low mo portion of the 
reach contour extends to m 1( / 2 ~ 1.3 TeV, and corresponds roughly to rrig ~ rrig ~ 3 TeV. 
The high mo portion of the reach contour extends to m\/i ~ 0.7 TeV, and corresponds to 
rrig ~ 1.8 TeV, while squarks are in the multi-TeV range, and essentially decoupled. Note 
that in this frame the parameter mo ranges all the way to 10 TeV. 

In frame b), with m^ = — 2.5mo, a much smaller range of mo is allowed, and the plot 
only extends to mo = 2 TeV. Since the reach is mainly determined by the values of the 
squark and gluino masses, we adapt the reach contours from Ref. |)l| to this non-mSUGRA 
case. Here, it is seen that the LHC reach covers almost all of the allowed A annihilation 
funnel. In frames c) and d), we show the cases for m^ = l.lm-o and 1.5mo, respectively. 
Here, the HB/FP type region re-emerges, but at much lower mo values, as discussed in 
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Sec. \Z.2[ The LHC reach is shown to cover all of the bulk and stau co-annihilation regions, 
but only a part of the higgsino annihilation (HB/FP) region. 




m D (TeV) m a (TeV) 



Figure 10: Approximate projections for the reach of the CERN LHC (100 fb _1 ) and ILC in 
the NUHM1 model in the mo vs. mi/ 2 plane for tan/3 = 10, Ao = 0, \i > and m t — 178 
GeV, for various choices of m^/m®. The regions shaded in green are consistent with the WMAP 
constraint flh 2 < 0.13, while those shaded in red and blue are respectively excluded by theoretical 
and experimental constraints discussed in the text. 

While we expect a similar reach of the LHC (in terms of rriq and rrig parameters) to be 
found in both the mSUGRA and NUHM1 models, the detailed gluino and squark cascade 
decays will change, as will the expected SUSY Higgs signals. To exemplify this, we list 
in Table [l] three model points. The first corresponds to mSUGRA for mo = wii/2 = 300 
GeV, Aq = 0, tan/3 = 10, fi > and mt = 178 GeV. We list a variety of sparticle and 
Higgs boson masses, along with 0.^ h , BF(b — > 57) and Aa M . The second and third points 
listed, NUHMla and NUHMlb, correspond to the same mSUGRA model parameters, but 



- 22 - 



with 77i0 = —735 GeV and 550 GeV, respectively. The mSUGRA case can be seen to have 
Q^/i 2 = 1.2, and is thus strongly excluded by WMAP data, while the two NUHM1 points 
have Q^h 2 ~ 0.11, and give the correct amount of CDM in the universe. The NUHMla 
point has a similar spectrum of sparticles compared to the mSUGRA case, although the 
heavier chargino and neutralinos have increased masses due to the larger value of the /i 
parameter. The main difference is that the heavier Higgs bosons are relatively light in 
the NUHMla case, and can be accessible to LHC searches as well as at a TeV-scale linear 
collider. In the case of mSUGRA, only the lightest Higgs h will be detectable at the LHC 
via direct h production followed by h — > 77 decay, or via tth or Wh production, followed 
by h — > bb decay. The h should also be observable in the sparticle cascade decays Q. In 
the NUHMla case with 777,4 = 265 GeV, the H and A Higgs bosons are much lighter, and 
should be detectable via direct H and A production followed by H, A — > rf decay p2| . 
The reaction gb — ► tH + followed by H + — > t + v t appears to be on the edge of observability. 
If tan (3 is increased to values beyond 15, then H, A —> fi + [i~ [|4|] should become visible. 



parameter 


mSUGRA 


NUHMla 


NUHMlb 


7770 


300 


-735 


550 


A* 


409.2 


754.0 


180.6 


m~ g 


732.9 


736.2 


732.0 


ma L 


720.9 


720.5 


722.4 


m h 


523.4 


632.4 


481.0 


m h 


650.0 


691.6 


631.0 


m iL 


364.7 


366.4 


364.5 


m~ eR 


322.8 


322.1 


323.0 


W2 


432.9 


759.6 


280.3 


m m 


223.9 


236.2 


150.2 


m z 4 


433.7 


759.5 


283.4 


m z a 


414.8 


752.0 


190.3 


m z 2 


223.7 


235.8 


160.7 


7775 


117.0 


118.7 


102.7 


m A 


538.6 


265.0 


603.8 


m H + 


548.0 


278.2 


613.0 


m h 


115.7 


116.1 


115.3 




1.2 


0.12 


0.11 


BF(b -» 87) 


3.2 x 10~ 4 


4.7 x 10~ 4 


2.5 x 10~ 4 




12.1 x 10~ 10 


9.4 x 1Q- 10 


17.4 x 10- 10 



Table 1: Masses and parameters in GeV units for mSUGRA and two NUHM1 models, where 
rno = my 2 = 300 GeV, A Q = 0, tan/3 = 10 and m t = 178 GeV. 

In the case of NUHMlb, is taken large enough that /i becomes small, 180.6 GeV, 
and the lightest neutralino develops a sufficient higgsino component to respect the WMAP 
dark matter constraint. The low /j value pulls the various heavier chargino and neutralino 
masses to low values ranging from 190 — 283 GeV. In this case, gluinos and squarks will be 
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copiously produced at the LHC. Gluinos dominantly decay via two body modes to t\t (BF ~ 
49%) and b 1>2 b (BF ~ 39%). The lighter stop decays via t x -► bW l)2 (53%), h -> tZ 3 (25%), 
while 61(62) mainly decays to the two charginos (roughly democratically to all charginos and 
neutralinos). The left squarks decay mainly to both charginos and to Z12, while q~R mainly 
decays via qn — * Z\ >2 . The lighter chargino decays via three body decays with branching 
fractions corresponding to those of the virtual W. On the other hand, Z 2 ^ decays via 
three body decays with the leptonic branching fraction BF(Z 2 ^ — ► iiZ\) ~ 1.5(3)% per 
lepton family. Finally the heavier chargino mainly decays via W 2 — > WZ 2 ,WiZ, while 
Z4 —* WiW. It is clear that the LHC will be awash in SUSY events, with gluino and 
squark production being the dominant production mechanism. In this scenario, the total 
SUSY cross section is almost 10 4 fb, so that even at the low luminosity, we should expect 
~ 100, 000 SUSY events annually. Moreover, from our discussion of the sparticle decay 
patterns, we see that all the charginos and neutralinos should be accessible via cascade 



decays of gluinos and squarks, as envisioned in Ref. |9J|. It would be extremely interesting 
to perform a detailed study of just how much information about the SUSY spectrum the 
LHC data would be able to provide in this case. While detailed simulation would be 
necessary before definitive statements can be made, it is plausible that analyses along the 
lines carried out in Ref. |)(| may yield information about a large part of the sparticle 
spectrum, and provide a real connection between collider experiments and dark matter 
searches. 

2.5.3 Linear e + e~ collider 

The reach of a = 0.5 and 1 TeV international linear e + e~ collider (ILC) for super- 
symmetry has been evaluated with special attention on the HB/FP region in Ref. [|97j in 
the case of the mSUGRA model. In this study it was shown that the reach contours in 
the mo vs. mi/ 2 plane are determined mainly via the reach for sleptons pairs, the reach 
for chargino pairs, and partly by the reach in Z\Z 2 production. There is also a significant 
reach for the Higgs bosons H, A and H + in the large tan (3 case. The striking result of 



Ref. [97] was that in the WMAP allowed HB/FP region, becomes small and charginos 



become light, the reach of the ILC extends beyond that of the LHC. In the HB/FP region 
of the mSUGRA model, squarks are in the multi-TeV regime, and effectively not produced 
at the LHC. The signal at the LHC becomes rate limited for nig ~ 1.8 TeV. 

In Fig. |l(], we also show contours of anc ^ m ^ n ( m e i ,e fl ) = 250 GeV and 500 GeV: 
since signals from chargino and selectron pairs can be probed at an e + e~ linear collider 
nearly up to the kinematic limit for their production, these contours follow the boundary of 
the region that would be probed at the ILC. Frame a) shows the mSUGRA case where an 
ILC would have an extended reach in the HB/FP region around mo ~ 8 — 10 TeV. In frame 
b) where the DM-allowed regions consist only of the A-funnel and stau-co-annihilation 
corridor, the ILC reach is well below that of the LHC. In frames c) and d), where > mo, 
the HB/FP region has moved to much lower values of mo- However, in these regions, again 
{l becomes small so that charginos become light, and the ILC reach for W^~W± pairs may 
exceed the reach of the LHC along the right-hand edge of parameter space, even in frame 
d) where squarks are comparatively light. 
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Concerning the specific SUSY models shown in Table |T|, in the case of the mSUGRA 
model, the ILC operating at y/s = 500 GeV would see of course Zh production, but also 
W^~W± and Z1Z2 production. The cross section for the latter process is ~ 200 fb, and 
since Z2 — > Z\Z® essentially all the time, the end points of the energy distribution of Z° 
should yield the values of Z\ and Z2 masses to good precision. This is, of course, over 
and above m^r r which can be determined as usual. In the NUHMla model, where now the 
MSSM Higgs sector becomes light, H°Z° and A°h production will also be possible, allowing 
a detailed study of the Higgs sector and possibly a good determination of tan f3 . If the 
ILC energy is increased somewhat above 500 GeV, then H + H~ also becomes accessible 
to study. For the NUHMlb model, the Higgs bosons again become heavy, but the various 
heavier charginos and neutralinos become light. In this case, the final states Z\Z%, Z\Z±, 
Z2Z2, Z2Z3, Z2Z/1 and even Z3Z4 as well as W^W^ are kinematically accessible. Thus, 
a whole host of heavier chargino/neutralino states would be available for study. If the 
ILC energy is increased to 1 TeV, then as in the mSUGRA case, the various slepton 
pair production as well as W1W2 pair production would be available for study and SUSY 
spectroscopy would become a reality. Moreover, the heavier Higgs bosons A, H and if* 
which now have substantial branching fractions to charginos and neutralinos, will also be 
accessible to study. 



3. NUHM2 model 



3.1 Overview 



The NUHM2 model is characterized by two additional parameters beyond the mSUGRA 
set. The two new parameters may be taken to be the GUT scale values of rn 2 H ^ and 
mjj , where these parameters may take on both positive and negative values. The model 
parameter space is given by 

mo, rn 2 Hu , m 2 Hd , m l/2 , A , tan/3, sign(^) . (3.1) 

We remind the reader that at tree level the Higgs scalar potential is completely specified 
by m H , ^ 2 an d the parameter Bfj,. The two minimization conditions allow us to 

trade two of nt^, m 2 H , fj. 2 and B\x in favor of tan(5 and M|, while a third may be traded 
for the CP odd Higgs scalar mass tua- In ( |3.1| ) above, fi 2 and B/j have been traded for 
tan/3 and M§, leaving the sign of \i (which enters via the chargino and neutralino mass 
matrices) undetermined. Alternatively, in the NUHM2 model, we could have eliminated 
m H u > m H d an d Bn leaving tan (3 together with the weak scale values of [i and as input 
parameters. Thus, the set 

m , /x, m A , mi/2, A , tan/?, (3.2) 

where fi, txia and tan/3 are input as weak scale values, while the remaining parameters 
are GUT scale values, provides an alternative parametrization of the NUHM2 model. In 
mSUGRA, we have two additional constraints, fn 2 H ^ = fn 2 Hd = rriQ, on the scalar potential 
and the values of fj, 2 and tua are determined. 
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We have upgraded Isajet v7.72 to allow not only the input of negative Higgs squared 
masses at the GUT scale, but also to accommodate the second of these parameter sets with 
weak scale values of fi and uia as inputs, using the non-universal SUGRA (NUSUG) input 
parameters p3[ . The NUHM1 and NUHM2 models in Isajet incorporate REWSB using the 
RG improved one-loop effective potential, minimized at an optimal scale Q = y/ m i L m i R 
to account for dominant two-loop contributions. 

An important aspect of the NUHM2 model is that RG running of soft masses is in 
general modified by the presence of a non-zero S term in Eq. (|2.8|) . The quantity S, which 
in the NUHM2 model is given by S = fn 2 H u ~ m H d ' en t ers the third generation soft scalar 
squared mass RGEs as, 

2 / 1 „ - n n 16 1 



I ^ g iMl-?, g lMl-^ g iMl + ^alS + f?X t + f b X b 



dt 16?r 2 V 15 3 10' 

(3.3) 



dm 



2 

~ b R — ^ I „2l\/r2 ^®„2l\/r2 i ^Je i o<2 



* - ifc. v u 9tMt - "3" + l 3ts+ Vi ' Xt )' (3 ' 5) 

The first and second generation soft mass RGEs are similar, but with negligible Yukawa 



coupling contributions. The Higgs boson soft mass RGEs are as given by (2.3) and ([2.4|). 
The coefficients of the S terms are all proportional to the weak hypercharge assignments, 
so that this term provides a source of intra-generational mass splitting. When S is large 
and positive (i.e., when m 2 Hu > m# ), the mass parameters for tr, en and jiR are the 
most suppressed, while those for qn and £l are enhanced. If S is large and negative, the 
situation is exactly reversed. For large values of \S\, the sfermion mass ordering as well as 
mixing patterns of third generation sfermions may be altered from mSUGRA expectations, 
or for that matter expectations in many other models of sparticle masses. For instance, 
it is possible that < ; moreover, while the lighter stau is usually expected to be 
dominantly tr in most models, this may no longer be the case in the NUHM2 model. 

As a simple illustration of how the spectrum of the NUHM2 model varies with the 
Higgs boson mass parameters, in Fig. O we show the physical masses of various sparticles 



versus Amu = uiq — sign(m 2 i [ ^) ■ w |?r^J = sign{m 2 Hd ) ■ y l^^/J — m o- This is a one 
parameter section of the NUHM2 parameter space we call the Higgs splitting (HS) model, 
with ArriH = corresponding to the mSUGRA model. Large positive Amn gives rise to 
a large negative S, and vice versa. In our example, we take uiq = mi/2 = 300 GeV, with 
Aq = 0, tan (3 = 10 and fi > 0. As Amn increases, we see that in the first generation, 
the en, da and ul masses all increase, while ei and ur masses decrease. In mSUGRA, 
rri£ R is always less than m^; in NUHM2 models, this mass ordering may be reversed. At 
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the highest allowed values of Arrin, the v T and f\ mass values become light (and fi is 
then dominantly enhancing i-channel Z\Z\ — > v T D T , e + e~ , VuVn, [i + fjT annihilation 
in the early universe, which lowers the relic density to within the WMAP bound. When 
mj> T , ~ rag , then co- annihilation reduces the relic density even further. 

m =300GeV, m 1/2 =300GeV, tan$=10, A =0, u>0, m t =178GeV 
;> 900 




-400 -200 200 400 600 800 1000 1200 1400 

Am H (GeV) 

Figure 11: Variation in sparticle masses versus Amu in the NUHM2 model, for mo = 300 GeV, 
m 1/2 = 300 GeV, A = 0, tan (3 = 10, fi > and m t = 178 GeV. 

Some aspects of NUHM2 phenomenology as a function of Amn are illustrated in 
Fig. [l^. In frame a) we show the values of //, ttia and mg versus Amu- For negative 
values of Aran, both fx and rriA are small, and we have a region of higgsino and (possibly) 
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Figure 12: Variation in a) tua, a* and , b) 0.% h 2 , c) BF(b — ► sj) and d) Aa^ versus Aran = 
(m Hti - m H „)/2 in the NUHM2 model, for m = 300 GeV, m 1/2 = 300 GeV, A = 0, tan/? = 10, 
/i > and tth = 178 GeV. 



A-funnel annihilation. For Aitih large and positive, the relic density, shown in frame 6), 
drops because left sleptons and sneutrinos become very light. The value of BF(b — > 57) 
and Aa^ are also shown in frames c) and c() . These rise for negative values of Aran because 
the charged Higgs bosons and the lighter charginos and neutralinos become very light with 
the -inos developing significant higgsino components. 

In Fig. [l^ we show again the variation in sparticle masses with ArriH, but this time for 
mo = 1450 GeV, m 1 / 2 =300GeV, Aq = 0, tan (3 = 10 and > 0. In this case, the large scalar 
masses yield a large S term in the RGEs, and the hypercharge enhancement /suppression 
is accentuated. We see, as noted in Ref. [41] for Yukawa unified models, that the ur and 
cr squarks are driven to very low mass values as Amu increases. At the high end of the 
ArriH range, they become the lightest squarks, even lighter than the t\. The large Arna 
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parameter space ends when ur becomes a charged/colored LSP, in violation of restrictions 
forbidding such cosmological relics. We also see that as Am,// increases, the t\ mass at first 
increases, then decreases, then increases again. The initial increase is because as Am// 
increases, X t decreases, leading to reduced Yukawa coupling suppression of the soft SUSY 
breaking mass parameters of top squarks. For still larger values of Am//, the S term grows 
and leads to a suppression of the t/j, and hence t\ mass. Finally, as Am// is increased even 
more, the Xt term becomes large and negative, and again resulting in an increase in the 
top squark soft masses. We have checked that throughout this range of Am//, the lighter 
top squark remains predominantly right-handed. 



In Fig. 14, we show the same frames as in Fig. |12|, but now for the mo = 1450 GeV case. 
At low Am// values, again fi gets to be small, so that the neutralino becomes higgsino-like 
leading to efficient annihilation in the early universe. At very large Am// values, the relic 
density is again in accord with the WMAP analysis - this time because the squarks become 
so light that neutralinos can efficiently annihilate via Z\Z\ — » uu, cc processes occurring 
via i-channel ur and cr exchange. Furthermore, if mn R and are in the 100-200 GeV 
range, they may be accessible to Tevatron searches! The light squarks also lead to a greatly 
enhanced neutralino-proton scattering rate, and hence to large rates for direct detection of 
relic neutralinos [JO]. There is a small enhancement in Aa^ at low Amu where charginos 
and neutralinos become light and higgsino-like, leading to larger chargino-sneutrino and 
neutralino-smuon loop contributions in the evaluation of (g — 2)^ [^] . 

3.2 NUHM2 model: parameter space 

Our first display of the parameter space of the NUHM2 model is in Fig. where we show 
the sign(m? H ^) ■ y'Wf/J vs. sign(m 2 Hd ■ yJ\fn 2 H \ plane for mo = m^j = 300 GeV, with 
Aq = 0, tan (3 = 10, /i > and rat = 178 GeV. In frame a) we show the allowed parameter 
space as the white region, while theoretically excluded parameter choices are red. The 
region to the right is excluded because /x 2 < 0. In the red region at the bottom, m\ < 0, 
while in that on the top, Z\ is not the LSP. The blue region is allowed theoretically, 
but here, dives to small values yielding < 103.5, in violation of the bound from 
LEP2. The parameter space of the NUHM1 model is shown by the black dashed line, where 
m # u = m ?f d ' wnne the mSUGRA value point where m 2 ^ = mjj = ra\ is shown by a black 
cross. The reader will notice that while the bulk of the parameter space of the NUHM2 
model lies above this dashed black line where "Z// d > m //„> there is a small portion for 
small values of \m 2 H | values this is not the case. The reason for this asymmetry is that 



as seen from the EWSB conditions ( [2, 10 ) and ( p. 11 ), the weak scale values of the Higgs 



mass squared parameters must satisfy fnj{ d > m H u > with m H u < : ^ the former inequality 
is badly violated at the GUT scale, radiative corrections cannot "correct this", and the 
correct pattern of EWSB is not obtained. We also show contours of [i (magenta) ranging 
from 300-2000 GeV, where 300 GeV contour is on the far right-hand side. Contours of rriA 
ranging from 300-1500 GeV are also shown, increasing from bottom to top. The -y/x 2 " value 
is shown in frame b), which shows most of the parameter space is excluded. The exception 
is the narrow green/yellow region near the lower edge of allowed parameter space, which is 
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SPS2: m^1450GeV, m 1/2 =300GeV, tan$=10, A =0, u>0, m t =178GeV 




1000 2000 3000 4000 5000 6000 



Am H (GeV) 

Figure 13: Variation of sparticle masses versus Am# defined in the text for the NUHM2 model, 
with to = 1450 GeV, m 1/2 = 300 GeV, A = 0, tan/3 = 10, [i > and m t = 178 GeV. 



the ^4-funnel, and on the right-most edge of parameter space, barely visible, is the higgsino 
region. The narrow green region at the upper boundary of parameter space corresponds 
to the slepton (or squark for large mo) co-annihilation region. 

The DM allowed regions of parameter space show up more prominently if the param- 



eters TTijj^ and Tnjj d are traded for tua and [i as inputs. We display in Fig. 16 the y x 2 
values for the same parameter space as in Fig. 15, but this time in the vs. [x plane. 
Again most of the parameter space is excluded, although in this mapping the higgsino re- 
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SPS2: m =1450GeV, m 1/2 =300GeV, tan$=10, A Q -0, \\>0, m t =178GeV 
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Figure 14: Variation in a) mA, M an d m z 1 i ^) ^z^ 2 ' c ) BF(b — * S7) and d) Aa^ versus Amj = 
(m Hti - mflJ/2 in the NUHM2 model, for m = 1450 GeV, m 1/2 = 300 GeV, A Q = 0, tan/3 = 10, 
/i > and rrn = 178 GeV. 



gion shows up as the broad band of green/yellow at low \i values, while the ^-annihilation 
funnel shows up as the vertical band running upwards near mA ~ 250 GeV. This plot 
highlights the importance of the measure of parameter space when deciding the likelihood 
that any particular framework satisfies some empirical (or, for that matter, theoretical) 
criteria: the tiny green/yellow sliver along the right edge in Fig. is expanded into the 
band, while the thin sliver at the bottom shows up as the A funnel. Notice also the thin 
green/yellow region for very large tjia values, where the existence of light sleptons brings 
the relic density prediction into accord with the WMAP value of VtcDhih 2 . The NUHM1 
model extends along the black dashed contour arc, with the mSUGRA model denoted by 
a cross. The regions away from the NUHM1 model arc denote SUSY mass spectra which 
are phenomenologically different from either the mSUGRA or NUHM1 model. In fact the 
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m =300 GeV, m t2 =300GeV, tar$=10, A Q =0, y\ >0, m, =178GeV 
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Figure 15: Plot of allowed parameter space in the vs. mn d plane of the NUHM2 model for 
for mo = m-i/2 = 300 GeV, Aq = 0, tan/3 = 10, fi > and mt = 178 GeV. In frame a), we show 
contours of /i and tua, while in frame b) we show values of x 2 . The yellow region in frame a) is where 
so-called "GUT stability bound" is violated. The short-dashed black line denotes the parameter 
space of the NUHM1 model. The cross denotes the mSUGRA model, while the long-dashed line 
gives the model where the Higgs scalar mass parameters are split as in Fig. [fl]. The blue region is 
excluded by LEP2. 

points with low /x and low itia, which only occur in the NUHM2 model, are somewhat 
favored by the combined constraints. The regions above and left of the GS contour violate 
the GS condition, while the HS contour denotes the path of the HS model through the 
NUHM2 model parameter space. 

The predictions for ^j^ft- 2 and contours for BFib — * s'y) and Aa^ are separately 
shown in Fig. [17] for a) fi > and also for the disfavored value b) fj, < 0. As before, we 
take rriQ = m\/2 = 300 GeV, Aq = 0, tan/3 = 10 and mt = 178 GeV. In frame a), the 
bulk of the allowed parameter space is determined by the WMAP allowed region. Within 
this region, the values of BF(b — > sj) and Aa M determine the best fit, which turns out 
to be rriA ~ 300 GeV and fx ~ 130 GeV. In this region, charginos and neutralinos as well 
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NUHM2: m =300GeV, m 1/2 =300GeV, tanp>=10, A =0, m t =178GeV 




m A (GeV) 



Figure 16: Plot of regions of \p\^ in the /i vs. rriA plane for mo = m 1 / 2 = 300 GeV, Aq = 0, 
tan/3 = 10 and m t = 178 GeV for fi > 0. The line labeled HS denotes the NUMH2 model where 
just the Higgs mass parameters are split as in Fig. |ll|. The region to the left of the Ah contour is 
where Ah production is accessible to a ^fs — 500 GeV ILC, while the region to the right of the f% 
contour is accessible to ILC via stau pair searches. 

as MSSM Higgs bosons are all relatively light. In frame b), for /x < 0, it is seen that 
BF(b — ► S7) is ~ 4 x 10~ 4 in the low right-hand region. For smaller values of \fi\ and rriA, 
the value of BF(b — ► 57) only increases, pushing the x 2 to large values all over the WMAP 
allowed region. 

3.3 Dark matter detection: the NUHM2 model 

We show in Fig. |l^ and [ll| the reach contours for the various detection channels introduced 
in Sec. |2.3| , respectively for the Burkert Halo Model (Fig. |l8|) and for the Adiabatically 
Contracted N03 Halo Model (Fig. ^): parameter space points lying below, or to the left, of 
the reach contours will yield detectable signals via the corresponding searches. The dashed 
black lines mark the locus of points appropriate to the one-parameter NUHM1 model, 
where m 2 H ^ = ni 2 Hd at Q = Mqut- The red cross indicates the particular point given by 
the universal mSUGRA case. 

For the chosen values of mSUGRA parameters (which include gluino and squark masses 
up to several hundred GeV) , stage-3 direct detection experiments will probe the bulk of the 
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NUHM2: m =300GeV, m 1/2 =300GeV, tan$=10, A^O, m,=178GeV NUHM2: m =300GeV, m 1/2 =300GeV, tan$=10, A^O, m,=178GeV 




250 500 750 1000 1250 1500 1750 2000 250 500 750 1000 1250 1500 1750 2000 



m A (GeV) m A (GeV) 

Figure 17: Ranges of f^/i 2 together with contours of BF(b — > sj) and Aa^ in the fi vs. tjia 
plane for juq — m 1 / 2 = 300 GeV, Aq — 0, tan/3 = 10 and mt — 178 GeV. For very large values of 
niA, we have the stau co-annihilation region. In frame a), we show contours for /i < and in frame 
b) we show contours for \i > 0. 



/i — mA plane allowed by cosmology, independently of the halo model under consideration. 
The exception is the region with very large values of mA ~ 1-8 TeV (not shown in the 
figure) where the relic density is in accord with the WMAP observation because the sleptons 
become very light. A large portion of the WMAP allowed region will also be within reach 
of stage-2 detectors, particularly for low values of the fi parameter. The shape of the VR 
curves is again readily understood in terms of the interplay of the two effects that we 
discussed in detail in Sec. 2.4: the enhancement of the heavy CP-even Higgs exchange 
channel at low mA, and the increased higgsino fraction, at low values of the \i parameter. 

The expected flux of muons from neutralino annihilations in the Sun is particularly 
sensitive to the value of the \x parameter: if the latter is sufficiently low, it provides a large 
enough spin-dependent neutralino-proton scattering cross section and a large capture rate 
of neutralinos in the core of the Sun, hence giving a large enough signal at IceCube. 



A comparison of Figs. [L8| and 19 shows what we had alluded to in the last section: 
the Burkert halo profile yields conservative predictions for the prospects for detection of 
relic LSPs. Indirect detection rates essentially track the size of the mass-rescaled pair 
annihilation rate (av)/m% . The difference between the profiles is especially accentuated 
for gamma rays, where the five orders of magnitude enhancement mentioned in Sec. 2.3 
implies that the entire plane will be covered for the Adiabatically Contracted Halo Model. 
The difference between the models is also considerable for the various antimatter searches, 
where once again, the Burkert profile leads to the most conservative prediction. 



We show in Fig. |2C] (m = 300 GeV) and Fig. [2J (m = 1450 GeV) the neutralino- 
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Figure 18: The reach contours (i.e. the iso-VR=l lines) for various dark matter detection tech- 
niques in the (m^,/i) plane. Points lying to the bottom-left of the lines are within reach of the 
future experimental facilities, as described in Sec. 2^3. The black dashed line indicates the param- 
eter space of the NUHM1 model in this plane, while the red cross locates the mSUGRA universal 
case. 



proton spin- independent scattering cross section a~ as a function of Am#, together 
with the current experimental limit and the projected sensitivity of Stage-2 and Stage-3 
detectors. In both cases, a~ increases significantly for large negative values of Am#, in 

Zip 

the region where the neutralino gets a large higgsino fraction since becomes small (left 
panels). Stage-3 direct detection experiments will thoroughly probe the WMAP allowed 
region in the left part of the plots, and Stage-2 detectors will have access to a wide portion 
of it. As for the large positive end of the Amn range, a WMAP compatible relic abundance 
is achieved, in both cases, through sfermion co-annihilations. In the thq = 300 GeV case 
squarks are heavy, and a light slepton sector does not particularly facilitate direct detection 
of the neutralino, which is found to lie even beyond Stage-3 detectors. 8 In the mo = 
1450 GeV case, the co-annihilating partners participate in the neutralino-proton scattering 
through s-channel squark exchange diagrams. This results in resonant squark contributions 
to neutralino nucleon scattering which enormously enhance a~ when nig — ► rn^. The 

steepness of a~ as a function of Amg. in the uiq = 1450 GeV case, is further increased 

Zip 

by the occurrence of an accidental cancellation between the up-squark neutralino-quark 
amplitude and the charm-squark gluon-mediated amplitude at Aitih — 5530 GeV. 

8 Although we can envision that the neutralino electron scattering may be resonantly enhanced if ~ 
mg i , the energy transferred to the electron is very small because the electron is very light. In this respect, 
the situation is different for the case of scattering off a nucleon discussed below. 
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Figure 19: The same as in Fig. but for the Adiabatically contracted N03 Halo Model. The 
VR contour for gamma rays from the galactic center is not shown, since the whole parameter space 
range shown in the figure features a VR larger than 1 in that detection channel (see the discussion 
in the text). 



Fig. [2^ (mo = 300 GeV) and Fig. |23] (mo = 1450 GeV) show our results for the ex- 
pected muon flux from the sun for the supersymmetric cases that we have been examining. 
As for the neutralino-proton spin-independent scattering cross section discussed above, the 
increase in the higgsino fraction in the left part of the left panels dictates a larger neu- 
tralino pair annihilation cross section as well as a larger spin-dependent neutralino-proton 
scattering cross section. As a result, when the neutralino thermal relic abundance enters 
the WMAP region, the expected muon flux gets close to or above the IceCube projected 
sensitivity for a hard neutrino spectrum, and even exceed the Super-Kamiokande limit in 
the m-o = 1450 GeV case. The abrupt drop in the muon flux visible in the mo = 1450 
GeV figure, close to the LEP-2 chargino mass limit, is due to the gauge bosons thresh- 
olds. As regards the sfermion co-annihilation regions (right panels), light sleptons induce 
only a slightly larger neutralino pair annihilation cross section (Fig. |22| ), while the squark- 
mediated resonant increase in produces an extremely large muon flux in the mo = 1450 
GeV scenario (Fig. p3| ). 

The overall summary and comparison of direct and indirect dark matter detection 
methods considered here is shown in Fig. |24| (mo = 300 GeV) and Fig. [25| (mo = 1450 
GeV) for the conservative Bukhert profile. This conservative assumption is why, in both 
cases, the gamma ray flux from the Galactic Center is well below the projected GLAST 
sensitivity. Antimatter searches will instead produce positive results in the higgsino-like 
region at negative Am//. In particular, antideuteron searches on GAPS look extremely 
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Figure 20: The neutralino-proton spin- independent scattering cross section a~ as a function 

Zip 

of Arnjj (black solid line), at m = 300 GeV. The other input parameters are specified in the 

figure. The dotted red line indicates the current 90% CL exclusion limit on a~ delivered by 

t-h ZlP 
the CDMS experiment |59J , while the blue dashed lines and the green dot-dashed lines correspond 

to the projected 90%CL exclusion limits for Stage-2 and Stage-3 detectors, taking as benchmark 

experiments the next-generation CDMS-11 detector and the XENON- 1 ton facility. 




promising in the present setup, as they will entirely probe the WMAP-allowed region for 
both values of mo- In the mo = 1450 GeV scenario (Fig. ^5|), as for the muon flux, 
antimatter fluxes and gamma ray flux undergo a sudden drop motivated by the gauge 
boson threshold (m^r < m\y). Light sfermions help to marginally increase the antimatter 
rates (see right panels), but the resulting fluxes are found to lie, in all cases, below the 
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Figure 22: The flux of muons from the Sun generated by charged current interactions of neutrinos 
produced in the center of the celestial body by neutralino-pair annihilations, integrated above a 1 
GeV threshold (solid black line), as a function of Amj (black solid line), at m = 300 GeV, with 
other parameters as specified in the figure. The red dotted curves indicate the current 90% CL 
Super-Kamiokande limits on the muon flux from the Sun, while the dashed blue lines the expected 
90% CL IceCube sensitivity. We corrected for the energy threshold mismatch in the computation 
of the projected IceCube exclusion curve, and we used a hard neutrino spectrum in the left panel 
and a soft neutrino spectrum in the right panel (see the text for details). 
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Figure 23: Same as in Fig. |22|, but at mo = 1450 GeV. 



expected experimental sensitivity. 
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Figure 24: The Dark Matter detection Visibility Ratios (VR) for various detection strategies, 
as a function of the GUT-scale Higgs mass splitting Am#, at mo = 300 GeV. The other input 
parameters are specified in the figure. When the VR is above the "Visibility Threshold" (i.e. the 
VR=1 line) the corresponding model is within the reach of the corresponding dark matter search 
experiment. 
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Figure 25: Same as in Fig.g4[ but at m = 1450 GeV. 
3.4 NUHM2 model: Collider searches for SUSY 

In this section, we illustrate aspects of the NUHM2 model that lead to collider search 
possibilities not possible in the NUHM1 model. We illustrate our points by showing in Table 
| three cases where the NUHM2 model can be in accord with the WMAP measurement, 
but which give rise to unique collider phenomenology. All parameter points in this table 
have the same values of m.i/ 2 = 300 GeV, Aq = 0, tan/3 = 10, fj, > and m; = 178 
GeV. For the first point, labeled NUHM2a, we take m = 300 GeV, with fx = 220 GeV 



and rxiA = 140 GeV. It occurs in the lower-left region of the plot in Fig. |16| and gives 
a relic density ^^_h 2 = 0.10. It is characterized by both a low fi and a low tua value, 
unlike the NUHM1 model, which must have one or the other small, but not both, to be 
in accord with WMAP. This point yields light higgsinos and light A, H and if SUSY 
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Higgs bosons. The second point, NUHM2b, has mo = 300 GeV as well, but has input 
parameters m 2 Hd = (1651.7 GeV) 2 and m 2 H ^ = —(1051.7 GeV) 2 . It is characterized by 
relatively light left-handed sleptons and sneutrinos, due to the large S term in the RGEs. 
Finally, we show NUHM2c, which has m = 1450 GeV with m 2 Hd = (7047.3 GeV) 2 and 
m 2 Hu = —(4147.3 GeV) 2 . It is characterized by the presence of very light ur and cr 
squarks. Notice that the right sleptons are heavier than all the squarks, and that f\ ~ fx 
is the lightest of the charged sleptons 9 . 



parameter 


NUHM2a 


NUHM2b 


NUHM2c 


mo 


300 


300 


1450 




220 


933.2 


3443.7 


m A 


140 


1884.6 


7765.1 


o 

m H d 


-(506.4) 2 


(1651. 7) 2 


(7047.3) 2 


mir 

ti u, 


-(263.5) 2 


-(1051. 7) 2 


-(4147.3) 2 


m~g 


726.4 


739.4 


807.8 


mu L 


'von f* 


r~7 ,\C\ A 

74U.4 


1 TO A O 

1724.8 


m UR 


7 lo.o 




i ri c 

151.6 


m h 


491.1) 


661.9 


1802.9 


m ' bl 


con r\ 

629.0 


730.6 


1830.5 


m e L 




lou.y 


OOU. 1 


m e R 


909 A 


0410. 


9^1 P. 1 


m fl 


290.1 


149.3 


522.9 


m„ T 


368.3 


129.9 


513.1 


m w 2 


293.8 


937.1 


3428.8 


m m 


174.4 


236.0 


250.4 


m z 4 


296.4 


935.7 


3427.1 


m z. 


228.5 


931.1 


3426.5 


m z 2 


178.7 


236.1 


251.5 


m z 1 


108.9 


119.2 


122.0 


m H + 


162.1 


1898.6 


7816.6 




113.3 


116.5 


120.3 


n^h 2 


0.10 


0.17 


0.14 


BF(b -► a 7 ) 


4.2 x 10~ 4 


3.3 x 10~ 4 


3.5 x 10~ 4 




15.3 x 10~ 10 


13.0 x 10~ 10 


1.2 x 10~ 10 



Table 2: Masses and parameters in GeV units for three NUHM2 models, where mi/2 = 300 GeV, 
A$ = 0, tan/3 = 10 and m t = 178 GeV. Input parameters are shown as bold-faced. 



9 There might also appear to be a possibility of generating characteristic spectra with S large and positive 
at the GUT scale, where only £r are light, while gluinos, charginos and neutralinos as well as most squarks 
and left sleptons would be heavy, making the signal difficult to detect at the LHC. This case does not, 
however, seem to be possible because large positive S leads to m 2 A < and thus a breakdown in the 
REWSB mechanism before the en becomes light enough. 
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3.4.1 Fermilab Tevatron 



The point NUHM2a with low [i and low vtla values will be difficult to probe at the Fermilab 
Tevatron. In this case, the chargino and neutralino masses are still rather large, and 
Z<i — ► Z\e + e~ has a branching fraction of just 1.5%, so that trilepton signals will be difficult 
to detect above background. The tua and tan/3 values are such that this point lies just 
beyond a "hole" in parameter space where none of the Higgs bosons are accessible to the 
Tevatron, even at the 95%CL exclusion level f9|, fUHj. The point NUHM2b, with relatively 



light sleptons, will also be difficult to probe at the Tevatron since charginos are quite heavy 
and slepton pairs are difficult to detect for any mass choices pOlfl . The point NUHM2c has 



two flavors of relatively light squarks - ur and cr - but the squark-neutralino mass gap 
is rather small. The signal would be identical to the one searched for in Ref. jl02j] , except 
with essentially twice the expected cross section for any given squark and LSP masses since 
o~(crcr) — o~(urur)). It is possible that a dedicated squark search might be able to detect 
a signal in the dijet+ Et channel, where relatively low jet Et and Et values ~ 25 — 50 
GeV might be expected, owing to the small ur — Z\ mass gap. Alternatively, if it becomes 
possible to tag c-jets with significant efficiency, it may be possible to suppress backgrounds 
sufficiently to pull out the signal. The phenomenology of light ur and cr squarks for the 
Tevatron is discussed more completely in Ref. plf . 

3.4.2 CERN LHC 

Squarks and gluinos would be produced at large rates at the CERN LHC. Their cascade 
decays would, in general, lead to the production of Higgs bosons in SUSY events. In the 
NUHM2 model, since the Higgs sector is essentially arbitrary, the heavier Higgs bosons 
H, A and could decay into other sparticles, resulting in characteristic events at the 



LHC [102]. In the case of the NUHM2a scenario shown in the table, the low value of ji 
implies the entire spectrum of charginos and neutralinos will be quite light, and accessible 
via squark and gluino cascade decays. The cascade decays will be much more complex 
than in a typical mSUGRA scenario, but as for the NUHMlb point discussed earlier, 
potentially offer a rich possibility for extracting information via a variety of mass edges 
that would be theoretically present. The feasibility of actually doing so would necessitate 
detailed simulations beyond the scope of the present analysis. In addition, in this scenario, 
the heavier Higgs bosons are also quite light, and in fact the heavier inos are able to 
decay into A, H and with significant rates. The production of H and A followed by 
H, A — > t + t~ and perhaps also to /i + , should be detectable at the LHC. In addition, 
Wh or tih production with W and one of the tops decaying leptonically should also be 
detectable. 

The point NUHM2b is characterized by light sleptons. In this case, qi cascade decays 
will be lepton-rich, since Z2 — > e^e with a branching ratio ~ 10%, while W\ — > II^ + lv 
essentially 100% of the time. The decay Z<i — > vv is also allowed, but in this case the 
sneutrino decays invisibly. Gluinos nearly always decay to uru or crc. The squarks then 
decay via ur — > uZ\ and cr — > cZ\, so that gg production will give rise to ~ 4 jet + Et 
events. 
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The point NUHM2c, characterized by light ur and cr squarks, will allow for squark 
production at very large rates. In addition, gg production followed by g — > uru or crc 
will give rise to 2 — 4-jet + $t events, since the jets from ur — > uZ% decay will be rather 
soft owing to the small ur — Z\ mass gap. The cascade decay events in this case will be 
lepton-poor since Z2 decays mostly to Z\h or uru, crc final states, and qi is heavy so that 
charginos are not abundantly produced via their decays. 



3.4.3 Linear e + e collider 

Any scenario similar to that represented by point NUHM2a would be a bonanza for the 
ILC. In this case, a y/s = 0.5 — 0.6 TeV machine would be able to access both chargino and 
all four neutralino states as well as the heavy Higgs bosons H, A and H^. The difficulty 
would be in sorting out the large number of competing reactions, but here, variable center 
of mass energy and beam polarization would be a huge help. A complete reconstruction of 



chargino and neutralino mass matrices may be possible [104]. The low value of tha should 
serve to distinguish this case from a NUHMla-like scenario. 

In the case of NUHM2b, the fi and slepton states would be accessible to early 
searches, along with W^W^ production. Beam polarization would be a key ingredient in 
determining that the fi and ei are left-handed. Determination that t\ is dominantly tl 
and/or mg L <C rng R would already point to an unconventional scenario. 

The case of NUHM2c would allow urur and crcr production to occur at large rates 
at an ILC. Again, the beam polarization would easily determine the right-hand nature of 
these squarks, which would be a key measurement. It should also be possible to determine 



their masses [105], and if it is possible to tag c jets with reasonable efficiency, to also 
distinguish between squark flavors. 

We note here that in addition, in the NUHM2 model, the reach of an ILC may be 
far greater than the CERN LHC for super symmetry. The reason is that the LHC reach 
is mainly determined by the mo and m^2 parameters, which determine the overall squark 
and gluino mass scales. In contrast, the ILC reach for chargino pair production depends 
strongly on the \x parameter. Thus, the NUHM2 case where m-o and mj/j are large, while 
H is small may mean chargino pair production is accessible to an ILC while gluino and 



squark pair production is beyond LHC reach. The case is illustrated in Fig. 26, where 
we show the mo vs. plane for A§ = 0, tan/3 = 10, fi > mt = 178 GeV and a) 

fi = rriA = 500 GeV and b) fj, = tua = 300 GeV. The yellow and green regions are WMAP 
allowed, while the unshaded regions have Qg^h 2 bigger than the WMAP upper bound. The 
yellow bands just above the LEP excluded blue regions in both frames is where 2m^ ~ nih- 
The corresponding band in the left panel at m^ 2 — 0.6 TeV is the A funnel, while in the 
upper yellow/green regions in both panels the LSP has a significant higgsino content. 

The SUSY reach of the CERN LHC should be similar to the case of the mSUGRA 



model calculated in Ref. [91|, and as before, we show this result as an approximate depiction 
of the LHC reach for the case for the NUHM2 model. We also show the mass contour in 
a) for a 250 and 500 GeV chargino, accessible to a \/s = 0.5 or 1 TeV ILC machine. Here, 
the 1 TeV machine has a reach beyond the large mo reach of the LHC. In the case of frame 
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b), fi is so small that the W\ mass is almost always below ~ 330 GeV, and so the entire 
plane shown would be accessible to a 1 TeV ILC! 



NUHM2: tanf,=10, A„=0, m A =500GeV, \i=500GeV, m t =178 GeV NUHM2: tanf,=10, A„=0, m A =300GeV, ^=300GeV, m t =178 GeV 
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Figure 26: Approximate reach of the CERN LHC (100 fb 1 ) and ILC for supersymmetric matter 
in the NUHM2 model in the too vs. mi/ 2 plane for tan/3 = 10, Aq = 0, /i > and m t = 178 GeV, 
for a) pi, = 500 GeV, m A = 500 GeV and 6) fj, = 300 GeV, m A = 300 GeV. 



4. Concluding Remarks 

We have examined the phenomenological implications of gravity-mediated SUSY breaking 
models with universal matter scalars, but with non- universal Higgs soft SUSY breaking 
masses. For simplicity, we assume a common GUT scale mass parameter for all matter 
scalars - this guarantees that phenomenological constraints from flavor physics are re- 
spected - but unlike in mSUGRA, entertain the possibility that the soft SUSY breaking 
mass parameters in the Higgs sector are unrelated to the matter scalar mass. In these 
non-universal Higgs mass (NUHM) models where the Higgs fields H u and Hd originate in 
a common multiplet (as, for instance, in an SO(10) model with a single Higgs field), we 
would have rn 2 Hu = m 2 H , and there would be just one additional new parameter (NUHM1 
model) |2l| vis a vis mSUGRA, while the more general scenario would have two addi- 
tional parameters (NUHM2 model). We have found that, once WMAP constraints are 
incorporated into the analysis, this seemingly innocuous extension of the mSUGRA pa- 
rameter space, which naively would not be expected to affect squark, gluino and slepton 
masses, significantly expands the possibilities for LHC and linear collider phenomenology 
from mSUGRA expectations: phenomena that were considered unlikely because they were 
expected to occur only in particular corners of parameter space become mainstream in the 
extended model. In the absence of any compelling theory of sparticle masses, the necessity 
for ensuring that all experimental possibilities are covered is sufficient reason to examine 
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the consequences of NUHM models, regardless of whether or not one considers these to be 
theoretically attractive. 

We have examined the allowed parameter space of the NUHM models with respect to 
neutralino relic density (WMAP), BF(b — > sj) and Aa M constraints. The WMAP upper 
bound on QcDMh 2 requires rather efficient LSP annihilation, and so severely restricts any 
supersymmetric model. Aside from the bulk region with small values of bino and scalar 
masses, enhanced LSP annihilation may occur if the LSP has significant higgsino or wino 
components (the latter is not possible in models with unified gaugino masses 10 ), resonantly 
annihilates via Higgs scalars (or Z bosons), or co-annihilates with the stau or some other 
charged sparticle. Within mSUGRA, the higgsino annihilation region occurs only if mo 
is very large, while resonance annihilation with heavy Higgs scalars is possible only for 
large values of tan (3. However, even in the simple one parameter NUHM1 extension of 
mSUGRA, for almost any values of mo, m 1( / 2 and tan/3, there are two different choices of 
irij) (defined in the text) that can bring the relic density to be in accord with the WMAP 
measurement: for large positive m^, one enters the higgsino region, while for large negative 
values of m^, one enters the A annihilation funnel. The higgsino region with small [i values 
gives rise to large rates for direct and indirect detection of neutralino dark matter, and also 
leads to light charginos and neutralinos which might be accessible to a TeV-scale linear 
e + e~ collider, or which can enrich the gluino and squark cascade decays expected at the 
CERN LHC. The A-funnel region in the NUHM1 model can also occur at any tan/3 value, 
and usually leads to relatively light H, A and Higgs bosons which may be accessible 
to collider searches. Also, the expected suppression of e and \x signals from cascade decays, 
which is expected in the mSUGRA model for points in the A-funnel due to enhanced -ino 
decays to taus[107|, will not necessarily obtain in the NUHM1 model since tan/3 is not 



required to be large. Since the early universe neutralino pair annihilation cross sections 
are enhanced on the ^-resonance, indirect DM signals are, in general, enhanced as well. 

The parameter freedom is enhanced even more in the NUHM2 model. In this case, the 
mSUGRA-fixed parameters \i and rriA can now be taken as inputs, rather than outputs. 
This allows one to always dial in a low value of fx such that one is in the higgsino region, 
or a low value of itia so that one is in the ^4-funnel. As before, direct and indirect DM 
detection rates are enhanced in these regions. Collider signals may change as well, since 
now all charginos and neutralinos can be light, and one can have enhanced cascade decays of 
squarks and gluinos to charginos, neutralinos and to heavy Higgs bosons. In the case where 
mo and m^j are large, but [i is small, the reach of a LC may exceed that of the CERN 
LHC. In the NUHM2 model, qualitatively new regions emerge where the relic density is 
suppressed due to novel sparticle mass patterns: very light left-handed sleptons, or very 
light right-handed up and charm squarks, which have obvious implications for collider 
signals. The latter case results in large rates for direct and indirect detection of neutralino 
DM, in addition to large jet+ fix signals at hadron colliders, and to the possibility of 
squark NLSPs at e + e~ LCs. 

In conclusion, we have seen that the seemingly innocuous decoupling of scalar Higgs 



For an exception to this, see the 200 model in Ref. JL06|. 
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mass parameters from other scalar masses can significantly alter our expectations of what 
we may expect in terms of dark matter as well as (s) particle physics phenomenology. This 
is partly because altering the Higgs potential can dramatically change the value of /j? 
that yields the correct value of M§, and partly because renormalization group evolution of 
sparticle mass parameters is dramatically altered by a non-zero value of S in the NUHM2 
model. For both the NUHM1 and NUHM2 extensions of the mSUGRA model one is able 
to find generic regions of parameter space that are in good agreement with the WMAP 
determination of the cold dark matter relic density, as well as with constraints from b — > sy 
and {g — 2) M . These regions can lead to distinctive signals at both direct and indirect 
dark matter detection experiments, and also provide distinctive signatures at both the 
CERN LHC pp collider and the International Linear Collider, with a center of mass energy 
^fs = 0.5 - 1 TeV. 
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